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Preface 


This  investigation  demonstrated  the  capability  of  the  Geometri- 
cal Theory  of  Diffraction  (GTD)  to  analyze  the  fields  due  to  mono- 
pole arrays  on  cylinders.  Although  one  of  the  principle  arguments 
for  the  use  of  the  GTD  method  to  solve  for  the  field  equations  was 
its  requirement  for  less  computer  usage  compared  to  alternative 
methods,  a computer,  as  such,  was  not  used. 

In  this  regard  I would  like  to  thank  Mr.  Robert  Barham  for  the 
use  of  his  HP-67  programable  calculator.  I would  also  like  to  thank 
Professor  Kouyoumjian  at  Ohio  State  for  his  comments  and  his  help 
in  obtaining  many  of  the  sources  on  GTD  used  herein.  I would 
especially  like  to  thank  Professor  Bill  Davis  for  taking  the  time 
to  go  through  most  of  the  theory  with  me  and  for  his  comments  on 
the  grammatical  problems  that  occured. 

Bruce  A.  Thieman 
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Abstract 


’Using  the  Geometrical  Theory  of  Diffraction  (GTD)  the  fields 
due  to  a monopole  array  mounted  on  an  aircraft  near  the  top  or 
bottom  were  analyzed.  The  aircraft  was  modeled  in  its  most  basic 
form;  an  infinitely  long  elliptical  cylinder.  The  aircraft  cross- 
section  at  the  location  of  the  antenna  was  modeled  as  a 2 meter  by 
4 meter  ellipse.  The  frequency  band  of  operation  was  250  mHz  to 
400  mHz  with  an  average  wavelength  assumed  of  1 meter.  The  three 
element  array  was  found  to  have  a beam  broadening  effect  when 
compared  to  the  array  mounted  on  an  infinite  ground  plane.  The 
array  exhibited  a poor  potential  to  form  a single  major  lobe  in 
its  antenna  pattern  but  had  a high  mobility  of  pattern  nulls.  .1 


ANALYSIS  OF  ANTENNA  ARRAYS  ON  CYLINDERS 


BY  THE 


GEOMETRICAL  THEORY  OF  DIFFRACTION 


I.  Introduction 


One  of  the  principal  problems  in  the  design  of  a reliable 
airborne  antenna  is  the  effect  on  volumetric  antenna  radiation  pat- 
terns by  the  aircraft  not  acting  as  a flat  ground  plane.  The 
particular  antenna  pattern  problem  to  be  considered  here  consists 
of  an  antenna  array  of  quarter-wavelength,  radial  monopoles,  which 
operate  in  the  250  mHz  "to  400  mHz  range,  mounted  on  an  aircraft 
fuselage  near  the  top  or  bottom.  The  array  configuration  is  depicted 
in  Fig.  1. 


Fig.  1.  Array  Configuration 


Since  this  is  a general  analysis,  the  aircraft  is  modeled  in 
its  most  basic  form.  The  fuselage  is  assumed  to  be  an  infinitely 
long  perfectly  conducting  elliptical  cylinder.  The  cylinder 


1 


solution  employed  in  this  analysis,  when  compared  to  the  experi- 
mental case  of  one  monopole  on  a cylinder,  gives  an  accuracy  of 
about +1  dB  in  pattern  comnutation  down  to  a level  about  20  d3 
below  the  pattern  maximum. 

Antenna  location  on  an  aircraft  is  usually  based  primarily 
upon  convenience  with  respect  to  the  aircraft  structure  with  no 
regard  for  system  requirements.  The  antenna  system  performance  is 
determined  by  analyzing  antenna  patterns  of  scale  models.  This 
approach  to  antenna  design  and  location  requires  large  amounts  of 
engineering  time  and  money.  This  investigation  develops  an  analytic 
approach  to  determining  the  antenna  system  performance. 

The  first  analytic  solutions  used  to  determine  on-aircraft 
antenna  patterns  were  the  cylindrical  harmonic  solutions  for  in- 
finitely long  circular  (Ref  1;  Ref  2)  and  elliptic  cylinders  (Ref  3). 
These  solutions  modeled  the  fuselage  by  a cylinder  whose  elliptical 
cross-section  approximated  the  fuselage  cross-section  at  the  antenna. 

'• 

location.  These  solutions  failed  for  structures  of  many  wavelengths 
in  diameter  unless  expanded  asymptotically  for  large  dimensions. 

The  expansion  of  these  solutions  required  computers  of  very  large 
memory. 

A more  recent  approach  is  the  use  of  the  integral  equation 
method  via  moment  methods.  By  solving  the  boundary  value  problem 
for  fields  on  aircraft  structures,  the  surface  currents  and  the 
resulting  scattered  fields  car.  be  found.  Richmond  (Ref  4;  Ref  5)> 
developed  a moment  solution  based  on  a wire  grid  modelling  structure. 
This  solution  required  the  determination  of  approximately  100 
unknown  currents  per  square  wavelength  in  order  that  the  grid 
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adequately  simulate  a perfectly  conducting  surface 


Another  approach  is  the  surface  current  model  method  (Ref  6). 
This  method  divides  the  surface  of  a conducting  body  into  patches, 
with  each  patch  having  two  orthogonal  unknown  complex  currents. 

This  reduces  the  number  of  unknown  currents  to  about  20  per  square 
wavelength  allowing  solutions  for  much  larger  surfaces.  This  solu- 
tion, modal  solution,  and  the  moment  method  are  restricted  to  lower 
frequencies  due  to  the  limitation  of  the  size  of  matrices  which 
modern  computers  can  solve  without  unacceptable  loss  of  accuracy. 

Another  approach  that  has  obtained  success  at  analyzing  on- 
aircraft  antennas  is  the  Geometrical  Theory  of  Diffraction  (referred 
to  simply  as  the  GTD  henceforth).  The  GTD  is  a high  frequency 
solution  which  is  divided  into  two  basic  problems;  wedge  diffraction 
and  curved  surface  diffraction.  Both  of  these  diffraction  solutions 
have  been  applied  in  computing  the  patterns  of  slot  antennas  mounted 
on  cylinders  (Ref  7;  Ref  8;  Ref9)i  rockets  (Ref  10),  and  aircraft 
(Ref  11;  Ref  12).  Using  this  approach,  one  applies  a ray  optics 
technique  to  determine  components  of  the  field  incident  on  the 
various  scatterers.  Components  of  the  diffracted  field  are  found 
using  the  GTD  solutions  in  terms  of  rays  which  are  summed  with  the 
geometrical  optics  terms  in  the  far  field. 

The  approach  aoplied  in  this  thesis  is  based  on  previous  work 
by  Burnside  (Ref  13)  and  Yu  (Ref  14),  which  demonstrated  the  capa- 
bility of  the  numerical  solutions  to  predict  the  radiation  patterns 
of  fuselage  mounted  antennas  in  an  efficient  and  economical  way.  If 
the  volumetric  patterns  were  found  directly  by  analyzing  rays  on 
complex  three-dimensional  surfaces  as  done  in  Reference  12,  the 
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resulting  numerical  solution  would  be  very  complex,  time-consuming, 
and  uneconomical.  Nevertheless,  if  certain  assumptions  can  be 
made,  the  approach  undertaken  previously  can  be  used  to  overcome 
these  difficulties  and  simplify  the  problem  a great  deal. 

The  basic  theoretical  background  on  the  geometrical  theory  of 
diffraction  (GTD)  for  curved  surface  diffraction  is  presented  in 
Chapter  II.  Expressions  for  the  parameters  and  special  functions 
involved  are  given  in  Chapter  III.  The  equation  for  diffraction  by 
a monopole  on  an  infinitely  long  cylinder  are  presented  in  Chapter 
III,  as  is  the  mathematical  model  of  the  aircraft.  In  Chapter  IV, 
the  numerical  results  obtained  for  the  single  monopole  case  are 
compared  with  published  results  and  measurements  and  in  addition 
the  field  patterns  of  the  three  element  array  for  different 
phases  are  presented.  Chapter  V contains  the  conclusions. 


k 
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II.  Theoretical  Background 


For  scattering  objects  that  are  large  in  terms  of  the  wavelength, 
such  as  for  electrically  large  aircraft,  the  GTD  is  the  optimum 
approach.  The  following  developement  of  the  GTD  assumes;  that  ray 
trajectories  are  determined  by  the  generalized  Fermat's  principle, 
that  the  field  amplitude  variation  along  a ray  path  can  be  found 
from  power  conservation  in  a tube  (or  strip)  cf  rays,  and  that  the 
diffraction  and  excitation  of  waves  can  be  treated  as  local 
phenomena.  The  validity  of  these  assumptions  and  the  GTD  have  been 
established  in  publish  d materials. 

According  to  geometrical  optics  the  region  exterior  to  a per^ 
fectly-conducting  surface  with  a radial  electric  current  source  at 
point  0 is  divided  into  an  illuminated  region  and  a shadow  region 
by  a plane  tangent  to  the  surface  at  called  the  shadow  boundary. 

A portion  of  these  two  regions  adjacent  to  the  shadow  boundary  is  a 
transition  rogion.  The  angular  extent  of  the  transition  region  is 
of  the  order  when  the  source  is  on  the  surface;  being 

the  radius  of  curvature  of  the  surface  at  Q in  the  direction  of 
propagation,  and  k being  the  wave  number  2 TT/A.  The  three  regions 
are  depicted  in  Fig.  2.  To  solve  for  the  fields  due  to  the  array, 
the  derivation  of  the  high-frequency,  far-zone  field  due  to  a point 
source  has  to  be  determined  for  these  three  regions* 


Diffraction  by  _a  Curved  Surface 

When  an  incident  ray  strikes  a smooth,  curved  perfectly  con— 
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ducting  surface  at  grazing  incidence,  i.e.,  at  the  shadow  boundary, 
a part  of  its  energy  is  diffracted  into  the  shadow  region.  This 
phenomena  is  described  by  a class  of  diffracted  rays  now  know  as 
creeping  waves  (Ref  15 )•  These  ray  paths  include  the  points  and 
Q2  which  form  a curve  on  the  diffracting  surface,  as  illustrated  in 
Pig.  3.  The  actual  concept  of  creeping  waves  was  introduced  by 
Franz  and  Depperman  (Ref  16;  Ref  17 )•  The  basic  concepts,  as  pre- 
sented in  the  following  discussion  is  taken  in  part  form  "Asymptotic 
High-Frequency  Methods"  by  Kouyoumjian  (Ref  18). 

The  diffraction  by  a smooth  curved  surface  is  shown  in  Fig.  3 
in  which  0 is  the  source  point  and  P is  the  observation  point  in  the 
shadow  region.  Diffracted  rays  are  characterized  by  an  extended 
form  of  Fermat's  principal  which  is:  The  diffracted  rays  connecting 
two  points  0 and  P in  the  exterior  of  a perfectly  conducting  body 
are  those  curves  joining  0 and  P which  have  stationary  (optical) 
length  among  all  curves  from  0 to  P having  an  arc  on  the  body  sur- 
face. Applying  this  principle,  the  line  OQ^Q^P  is  the  shortest 
distance  between  0 and  P which  does  not  penetrate  the  surface.  In 
detail,  a ray  incident  on  the  shadow  boundary  at  divides;  one 
part  of  the  incident  energy  continues  straight  on  as  predicted  by 
s-eometrical  optics,  a second  part  follows  the  surface  S into  the 
shadow  region  as  a surface  ray  shedding  diffracted  rays  tangentially 

A A A 

as  it  propagates  where  t,  n,  and  b are  the  unit  vectors  in  the 
direction  of  incidence,  normal  to  the  surface  S and  binormal  to  the 
surfave  (b=t  x n).  The  incident  field  E (Q^)  may  be  seperated  into 
its  normal  and  binormal  comronents  (n  * B1(Q^)  and  b • E^(Q^)). 

It  is  assumed  that  these  two  components  induce  surface  ray  fields 
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top  view 
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Fig.  3.  Diffraction  by  a smooth  convex  surface 
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which  propagate  independently  of  each  other  along  the  geodesic  arc 
between  and  Qg. 

The  binormal  surface  ray  field  at  is  related  to  the  binormal 
component  of  the  incident  field  at  by 


f 


I 


L. 


A(Q1)exp( jf(Q1))  - DS(Q1)(b1  • E1^))  (l) 

where  DS(Q^)  is  the  scalar  diffraction  coefficient  for  a soft 

boundary  condition.  The  terms  hard  and  soft  boundary  conditions 

are  terms  from  acoustics.  A soft  boundary  is  one  where  the  preasure 

field  is  zero  at  the  boundary  (E  = Ez?).  A hard  boundary  condition 

is  one  in  which  the  derivative  in  the  direction  normal  to  the 

boundary  is  zero  ()e/^o  = 0}  H = H ?)•. 

z 

The  amplitude  of  the  surface  ray  is  assumed  to  be  governed  by 
the  conservation  of  energy  between  a pair  of  adjacent  surface  rays. 
This  principal  states  that  in  time  harmonic  fields,  the  energy 
flux  is  the  same  through  every  cross  section  of  a tube  of  rays,  such 


S’ig.  4.  Astigmatic  Bundle  of  Rays 
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2 2 
Let  A be  the  intensity  of  the  field  at  S and  Aq,  the  intensity 


at  the  reference  point  0.  For  an  isotropic,  homogeneous  medium 


2 2 
A d<T  = A do- 


(2) 


where  do-  is  the  cross  section  of  the  tube  of  rays  at  S and  do^  is 


the  cross  section  at  0.  Therefore 


A = A 


,Vd<rc/d<: 


(3) 


If  s is  the  distance  from  0 to  S,  then 


do- 


P l/2 


do-  + s)(^  + s) 


(4) 


where  ^ and  are  the  principal  radii  of  curvature  of  the  wave- 
front  through  0. 

Applying  the  energy  conservation  principle  to  a narrow  strip 
of  rays  on  a curved  surface  (see  Fig.  5)  gives 


A 2(J  + dJf)d dJP)  - A2(^)drj(^)  = -20((^)A2(>e)d^(i)d! 


(5) 


where  2<(4)  is  a proportionality  factor  relating  the  change  of  flux 
to  the  energy  radiated.  Dividing  by  djjand  letting  d|  tend  to  zero 
gives 


— (A2d^)  = -2of(A2d^)  (6) 
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Fig.  5.  A Band  of  Diffracted  Rays 
on  the  Surface  of  a Diffracting  Object 


Upon  integrating  Eq.  6,  the  result  is 


A(Q2)  = A(Qi 


Q2  1 
- ^(j?)dje 


where  d^  and  d^  = the  separation  between  a pair  of  rays  at  and 

Qg  respectively. 

- the  attenuation  constant  which  is  a function  of 
X and  the  coordinates  along  the  surface  and  de- 
pends on  the  local  radius  of  curvature  and  its 
derivatives. 

From  Fig.  3*  Q2  is  a caustic  of  the  diffracted  field  and  the 
second  coustic  is  located  a distance  from  Q^.  The  caustic  of  a 
given  family  of  rays  is  defined  as  that  surface  or  curve  which  is 
everywhere  tangent  to  the  ray  surface.  Thus  the  surface  of  a 
diffracting  body  is  a caustic  of  the  family  of  diffracted  rays  re- 
sulting from  a given  field  being  incident  on  the  body.  Thus  the 
binormal  component  of  the  diffracted  field  which  radiates  from  Q0 

c_ 

towards  P can  be  found  using  a caustic  at  for  a reference  point  as 


b2  • Ea(p)  = DS(Q2)A(Q?)exp( jfc2) 


(/*  + s) 


jxp(-jks) 
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where  D^Q^)  is  a diffraction  coefficient  of  the  soft  boundary. 
Based  on  the  reciprocity  principle — that  a source  at  0 produces 
the  same  field  at  P as  a source  at  P produces  at  0 — D^Q^)  has  the 
same  form  as  DS(Q^).  Combining  Eqs.  (l),  (7)>  and  (8)  yields 


The  incident  field  b^  • E (Q^)  excites  an  infinity  of  surface 
ray  modes,  each  with  its  own  diffraction  coefficient  and  attenuation 
coefficient.  For  example,  in  the  canonical  problem  in  Appendix  A, 
if  the  field  scattered  from  a cylinder  given  by  Eq.  (58)  is  approx- 
imated asymptotically,  it  is  seen  that  there  is  a surface  ray  mode 
associated  with  each  term  of  the  series,  i.e.,  with  each  value  of  m. 
Therefore,  Eq.  (9)  is  replaced  by 


This  equation  relates  the  diffracted  field  at  P to  the  incident 
field  at  for  the  soft  boundary  condition. 

An  expression  similar  to  Eq.  (10)  is  also  obtained  for  the 
n-rmal  component  of  the  incident  field: 


12 


with  the  superscripts  h and  s,  respectively  added  to  and  o^. 

The  diJ^j'd^  and^Oare  found  by  differential  geometry  involving  the 
rays  and  surfaces  this  is  discussed  at  length  in  Levy  and  Keller 
(Hef  15). 

To  obtain  the  diffraction  and  attenuation  constants,  the 

exact  canonical  solution  of  the  boundary  value  problem  for  the 

diffraction  in  question  must  be  known.  Such  a solution  is  expanded 

asymptotically  for  ka  large  where  a is  the  radius  of  curvature 

for  cylindrical  diffraction.  By  comparing  this  expansion  with  the 

GTD  solution,  the  values  of  ot  and  D can  be  determined.  Appendix  A 

mm 

contains  a derivation  of  the  cerer'-.lizec  diffraction  coefficient 


" ^ 

and  attenuation  coefficient  for  diffraction  by  curved  surfaces. 

The  diffraction  thus  far  discussed  is  applied  to  an  open 
curved  surface.  For  a closed  surface,  each  surface  ray  mode  pro- 
duced at  encircles  the  surface  an  infinite  number  of  times.  The 
length  of  the  surface  ray  path  for  the  nth  encirclement  is  + nT 
where  T is  the  distance  traveled  in  one  transversal  of  the  closed 
surface.  These  multiple  encircling  rays  may  be  summed  to  contribute 


to  the  denominator  of  the  diffracted  field.  Note  for  this  derivation 
that  the  exp(jkz)  convention  has  been  adopted. 


Direct  Rays 

The  rays  emanating  from  a point  current  element  with  current 
moment  Iax  determines  the  incident  field  for  the  diffracted  and 
reflected  rays  as  well  as  the  field  at  all  points  where  Fermat’s 
principle  produces  a straight  line  from  the  source. 

Referring  to  Fig.  6,  A is  defined  as  (Ref  19s  16 ) 

I A*  exp(-jks') 

Ax=-^ (15) 

4lTs' 

Applying  Maxwell's  equations 


H - V x Axx  - yAx  x £ U6) 


Pig.  6.  Geometry  for  current  source  at  origin 


Since 


then 


r X x * r X (rcosd  - S^inO)  = - ^sin©  (19 

H = + -^Jexp(-;jks')sineif'  (20 


2 

In  the  far  zone  field  the  l/s'  term  is  dominant,  the  1/s'  term 
tending  to  zero  much  faster  than  l/s'  therefore  it  can  be  neglected 
The  electric  field  is  obtained  from 


E 


7 X H 


(21 


xp(-jks 


sin0© 


+ 


to  give 


Again,  in  the  far  zone  field  the  l/s'  term  is  dominant  and  the 
direct  electric  field  becomes 

£ - ^jfexp(-jk8')  ^ sineS  (23) 

Launching  Coefficients  and  Constants 

The  theory  presented  up  to  now  has  assumed  a point  source  near 
a cylinder.  Monopoles  mounted  on  a cylinder  could  be  modeled  as  a 
line  of  point  sources  with  sinusoidally  distributed  amplitudes, 
using  this  theory.  The  far-zone  field  would  then  be  the  result 
of  vectorially  summing  the  fields  from  each  source  point  at  each 
field  point.  The  result  is  a. large  number  of  ray  tracings  and  field 
calculations  for  even  a simple  antenna  array. 

To  simplify  the  calculations,  the  sum  of  the  incident  fields 
due  to  all  of  the  point  sources  making  a monopole  model  are  approx- 
imated by  a launching  coefficient.  To  determine  the  launching 
coefficients,  the  radiation  from  an  axial  mounted  electric  current 
source  on  a perfectly  conducting  cylinder  is  analyzed.  The  asym- 
ptotic solution  of  this  canonical  problem  is  described  in  Appendix 
B,  where  it  is  found  that 

l/2  0 

Lm  = "("I")  H^2)'(ka)D®(Qi)  (23) 

m 

Lm  = -j(J|-)1/2H(.2)(ka;D^(ft1) 
m 

where  y are  the  zeroes  of  the  Hankel  function  in  the  first  equation 
m 

and  the  zeroes  of  *he  derivative  of  the  Hankel  function  in  the 
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i 

£• 


i 

i 


t 


second,  equation.  Note  that  is  now  the  source  location.  These 
launching  coefficients  are  valid  for  monopoles  whose  length  is  small 
in  comparison  to  the  dimensions  of  the  scatterers. 

Eq.  (12)  now  becomes 


E(P)  = (n2  v(Q1,Q2)  + b2  u(Q1,Q2))]C^2Ijexp(-jks)  (24) 


where  v(QltQ2)  and  u(Q^,Q2)  are  equal  to 


"*  oo  ' / Qg 

iexp(-jkl)^L  (Qi)D  (Q2)exp  - I (-f)d^  (25) 

2 m=0  m m ^ J o m J 


d*?2 


with  the  superscripts  h and  s,  respectively  added  to  Dm,  Lm,  and 


Diffraction  Coefficients  and  Attenuation  Constants 

The  diffraction  coefficients  and  attenuation  constants  depend 

on  the  local  geometry  of  the  surface,  the  wave  number  k,  and  the 

nature  of  the  surface,  as  described  by  the  boundary  conditions. 

Voltmer  ( Ref  20)  has  derived  the  diffraction  coefficients  and  attn- 

uation  constants  for  general  convex  surfaces  and  are  included  in 

Table  I.  Appendix  A contains  the  derivation  of  D and  al  for  a 

mm 

cylinder. 

To  first  order  these  functions  are 


exp(  .iTT/6 ) 


^5/V/5(^ji/6 


Pj- 


1/2 


'(-qm)); 


(26) 

(27) 
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CEHEMLUEO  DIFFRACTION  COEFFICIENTS  AND  ATTENUATION  CONSTANTS 


0 ■ radius  of  curvature  along  the  geodesic 

radius  of  curvature  perpendicular  to  tlto  geodetic  (traniveno  curve) 
Dots  Indicate  differentiation  vlth  respect  to  art  length  variable 


w 


for  the  "soft"  surface,  where  p is  the  radius  of  curvature  of  the 
surface  along  which  the  surface  ray  is  propagating  in  the  plane 
containing  the  normal  and  tangent  to  the  surface  ray , and 

kl/3 


exp(  -.iir/12 ) 

^ g 


1/2 


(D")  ‘^170=  Ui<^)): 


(28) 

(29) 


for  the  hard  surface,  where  the  Miller-type  Airy  function  is  given 

by  r 

Ai(-x)  » ~ J cos(-jt3  - xt  )dt  (30) 


Ai(-qm)  = 0 

Ai(“V  = 0 


(31) 

(32) 


where  q and  q are  the  zeroes  of  A.  and  A. ' respectively  with  the 
prime  denoting  differentiating  with  respect  to  the  argument  of  the 
function.  These  zeroes  are  calculated  in  Appendix  C. 

The  terms  "hard"  and  "soft"  refer  to  the  boundary  conditions. 
A soft  boundary  condition  is  a Direchlet  boundary  (the  field  van- 
ishes at  the  surface).  A hard  boundary  is  a Neumann  boundary 
condition  (the  normal  derivative  of  the  field  vanishes  at  the 
surface).  This  assumes  two  types  of  surface  ray  modes;  one  in 

— A 

which  E b bS,  satisfying  a soft  boundary  condition,  and  one  in 
which  E - nE,  satisfying  the  hard  boundary  condition.  The  cor- 
respondin'-.* coefficient  approximations  are  essentially  what  Keller 
and  Levy  (Ref  15:166-169)  derived.  Table  I contains  Voltmer's 
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corrected  results 


k 


L 


Transition  Regions 

The  series  representations  are  rapidly  convergent  when  the 
field  point  is  deep  in  the  shadow  region.  When  the  local  average 
radius  of  curvature  of  the  surface  is  larger  than  a wavelength  at 
the  points  of  diffraction,  only  the  1st  terms  are  neccessary  to 
achieve  good  accuracy.  As  the  field  diffraction  point  Qg  approaches 


(i.e.  Jf  becomes  very  small),  more  terms  must  be  added  to  maintain 


accuracy  and  it  becomes  inconvenient  to  treat  the  excitation, 
propagation  and  diffraction  of  the  surface  ray  seperately.  The 
series  representation  can  now  be  approximated  by  integral  represent- 
ations, which  have  been  found  to  be  proportional  to  Fock  type 
functions.  When  the  source  is  on  the  surface,  the  angular  extent 


of  the  transition  region  from  the  shadow  boundary  is  nearly 

1-1/3 


(kA>  J'/ J radians.  The  expressions  for  the  fields  in  the  transition 
region  are  also  deduced  from  the  asymptotic  solutions  to  appropriate 
canonical  problems  as  described  in  Appendix  E. 

The  well  tabulated  Fock  functions  (Ref  21)  are 
l)  Hard  Fock  function 


«((•  ) 


I | exp(-,j‘i 

w'(T 


zinl 

) 


dr 


(33) 


2)  Soft  Fock  function 


-(<;)  = ±,\ 


w2(T) 


(34) 


J 


which  are  described  in  Appendix  E and  tabulated  in  Appendix  G. 
The  argument  of  the  Fock  function  to  the  first  order  is 


M'2' 


When  the  field  point  is  in  the  shadow  region,  J > 0 and  when  the 
field  point  is  in  the  illuminated  portion  of  the  transition  region, 

X < 0,  where  X is  the  distance  traveled  on  the  surface  of  the 
cylinder.  Note  that  for  these  expressions  that  the  surface  ray 
divergence  factor  d^/d*^  is  equal  to  one  in  this  region. 

In  the  transition  region,  for  sources  very  close  to  the  surface 


the  term 


(Q2 

Iv-H 


in  Eq.  (25)  is  replaced  for  v(QltQ2)  by 


rg(0,  ^<0 

X>° 


and  for  u(Q1,Q2)  by 


where 


• s/5)g($),  t < 0 

J >0 

■ - f^)1/3 

i - </V«2>//^1>>1/6 
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III.  Diffraction  by  an 
Infinitely  Long  Elliptic  Cylinder 

The  special  case  of  the  GTD  solution  treated  here  is  one  in 
which  l/4A  monopole  antennas  are  mounted  on  the  curved  surface  of 
an  infinitely  long,  2A  by  4A  elliptical  cylinder.  An  average 
wavelength  of  1 meter  is  assumed  for  this  problem. 

The  particular  array  to  be  considered  is  depicted  in  Fig.  1. 
To  solve  the  field  problem,  each  monopole  is  approximated  by  a 
launch  coefficient.  The  field  at  the  field  point  ? will  be  the 
sum  of  the  fields  produced  by  each  source. 

Ray  Approximations 

Previous  work  by  Pathak  and  Luebbers  (Ref  22)  has  shown  that 
in  the  illuminated  region  (P^,  see  Fig.  7)>  the  diffracted  rays 
and  s,-,  have  a negligable  amplitude  compared  to  the  direct  ray  s^. 
For  point  sources  on  the  cylinder,  the  field  can  be  modeled  as 
only  ray  s^  and  Eq.  (23). 

In  the  shadow  region  only  the  diffracted  rays  s^  and  s^ 
contribute  to  the  field  at  P^.  Since  the  rays  leaving  (or  G) 
and  encircling  the  cylinder  n times,  cover  a cylinder  circumference 
of  about  nine  meters,  they  are  heavily  attenuated  and  their  effect 
is  very  small.  The  field  at  is  then  the  sum  of  the  rays  QjQ-2^3 
and  applied  to  Eq.  (24). 

In  the  transition  region,  the  integral  form  in  Eq.  (36)  and 
(37)  are  used  to  solve  for  the  field  at  F~ , avain  simplified  by 
ignoring  those  rays  which  encircle  the  cylinder  or.e  or  more  times. 
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Rays  Emanating  Fa0m  ^ Source  On  A Curveo  Surface 


Note  that  for  s^  the  diffraction  field  formula  in  Eq.  (23)  is  used. 
For  the  equations  in  all  of  these  regions,  the  spread  factor 
dijj/di^,  is  one. 

Field  Equations 

The  GTD  solutions  for  the  radiated  electric  field  of  a point 
source  on  an  elliptical  cylinder,  as  shown  in  Fig.  8,  with  all 
amplitudes  combined  into  A and  neglecting  torsional  effects  are 
Lit  Region 


E ■ -A  8in(6m)exp(-j(ks'  + <f)  0 


(40) 


Transition  Region 
a)  Lit  side 


E - |sin(©m)exp[-j(k(s  -J)  +«^g*(-$) 


(41) 


b)  Shadow  side 


E * f exp[-j(k(s  +i)  + (f)Jg *fc)^t(  " 


(42) 


Deep  Shadow  Region 


E = ^.A  exp£-j(k(s  + i)  + <f)J  n\ 


j 


T~  exp 
m=l 


(43) 
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Fig  8.  Geometry  of  nonopole  antenna  mounted  on  an  infinitely 
long  elliptic  cylinder 


Where 


complex  conjugate  of  the  Fock  function  (Ref  21) 

normal,  binormal  and  tangent  unit  vectors  to 
the  surface 

longitudinal  radius  of  curvature 

defined  in  Eq.  (35) 
phase  of  source 
amplitude  of  source 

distance  from  source  to  diffraction  point  or 
field  point  if  cylinder  is  not  intersected, 
distance  from  diffraction  point  ^ to  field 
point 

Geometry  of  Elliptical  Cylinder  ( Equations ) 

In  order  to  analyze  this  problem  by  GTD,  an  efficient  solution 
for  the  geodesic  paths  on  the  elliptic  cylinder  must  be  found.  The 
coordinate  system  adopted  here  is  illustrated  in  Fig.  9 and  defined 
b,r 

x = d cosh  u cos  v = a^.  cos  v 
y * d sinh  u sin  v * b^,  sin  v 
z = z (44) 

where  2d  is  the  distance  between  the  foci  of  the  ellipse.  To  define 
ar.  ellipse,  u must  equal  tanh  ^(b^/a^).  The  coordinate  v then 
defines  any  point  on  the  ellipse  and  varies  from  0 to  2tt. 

rsing  the  calculus  of  variations,  the  geodesic  paths  on  an 
ellirtic  surface  a-a  given  by 


«*(•) 


A ^ A 

n,  b,  t 


fc 


A 
s 1 
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Fig.  9»  Elliptic  Cylinder  Coordinate  System  (Ref  14:26) 


If  one  defines  the  geodesic  starting  direction  by  the  angle  (' 

as  shown  in  Fig.  8,  then  C = -cos  ei  . The  Parameters  of  the 

s 

problem  are 


z 


-cos4. 

s 

| sin*  | 


, o .2 

sm 


v ♦ b~,  cos  v dv 


(46) 
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1 

sin* 


2.2  ,2  2 , 

|af  sm  v + b^  cos  v dv 


-a^.  sin  vx  + b^cos  vy 

V~2  ~~i  ~2  r 

a ^ sin  v + b.p  cos  v 


(arc  length) 


(curvalinear 

coordinates) 


£ ^ A . 

t = sm*  e,  — cos*  e- 
si  s 2 


(unit  tangent 
vector) 


, A A 

^ b^.  cos  vx  + a^  sm  vy 

n = - ■ ; t 

•i/l  ~2 72  2 

|faf  sm  v + b^.  cos  v 


(unit  normal 
vector) 


A A ^ , A . A 

b a t X n = -cos*  e.  - sm*  e0  / . . 

si  s 2 (unit  bmormal  (52) 


vector) 


,2.2  2 2 n3/2 

(af  sm  v + bf  cos  v)  ' 


, • 2 
am-  sm  o< 


( longitudinal 
radius  of  curvature) 


Using  the  preceeding  Eqs.  (40)-(53)»  the  total  radiated  field  for  a 
ooint  source  on  a cylinder  can  be  determined. 


r 


K" 

Sit'- 
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IV.  Pattern  Calculations 


The  field  equations  (40)  - (53)  are  used  in  the  following 
figures  to  calculate  the  field  due  to  the  array  in  Fig.  1.  The 
tabulated  fields  due  to  each  monopole  are  included  in  Appendix  G. 

To  verify  the  general  GTD  equations  derived  for  a monopole, 
the  antenna  pattern  for  a single  monopole  on  the  cylinder  was  cal- 
culated and  compared  to  experimental  results  obtained  from  published 
literature  (Ref  3:662).  Fig.  10  contains  a plot  of  both  the  experi- 
mental data  and  the  calculated  values.  The  variation  of  the  experi- 
mental field  is  do  to  interference  from  the  rays  that  encircle  the 
cylinder  more  than  once.  It  was  assumed  in  the  calculations  that  for 
fields  not  in  the  deep  shadow  region,  that  only  the  dominant  term 
would  be  considered. 

Fig.  11  and  Fig.  12  (Ref  14:64)  provide  a means  of  comparison 
of  the  calculated  values  with  those  obtained  from  a scale  model  test 
of  an  aircraft.  Fig.  10  and  Fig.  11  are  almost  identical  except 
that  the  theoretical  case  exhibits  a stronger  field  in  the  deep 
shadow  region.  The  effect  of  wings  (Fig.  12)  on  the  aircraft  causes 
an  even  further  degradation  of  field  strength  in  the  shadow  region. 
The  calculated  antenna  patterns  are  therefore  valid  as  a good  approx- 
imation for  0 within  the  range  of  + 105°  with  respect  to  the  antenna 
when  the  antenna  is  mounted  near  the  top  or  bottom  of  the  fuselage. 

Fig.  13  is  a plot  of  the  antenna  pattern  of  antenna  2 or  3.  Do 
to  symmetry,  their  patterns  are  reflections  of  each  other.  The 
effect  of  mounting  the  antenna  off  of  the  main  axis  of  the  cylinder 
is  to  shift  the  nulls  of  the  pattern  as  well  as  cause  sharper 
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changes  in  the  pattern  where  the  minima  of  the  radius  of  curvature 
of  the  cylinder  occur.  Note  that  the  null  has  shifted  about  26° 
from  its  position  in  Fig.  10  although  antenna  2 and  3 are  mounted 
perpendicular  to  the  cylinder  surface  at  a point  that  is  about  14° 
off  of  the  primary  axis  of  the  ellipse  (v^  «■*  14°).  This  shift  is 
caused  by  the  antenna  angle  in  space  being  about  26°  from  the  ver- 
ticle. 

Fig.  14  is  a plot  of  the  antenna  pattern  with  all  three 
elements  of  the  array  driven  with  the  same  amplitude  and  phase. 
Comparing  this  pattern  to  that  of  a 3-element  triangular  array  de- 
picted in  Fig.  15*  driven  under  the  same  conditions,  shows  that  the 
cylinder  has  the  effect  of  broadening  the  lobes  as  well  as  increasing 
the  angle  • where  the  pattern  maxima  occurs.  The  result  is  a pattern 
more  like  that  of  a single  monopole  or  dipole  in  space. 

Fig.  16  is  a plot  of  the  antenna  array  pattern  with  element  1 
driven  45°  out  of  phase  with  elements  2 and  3.  It  takes  very  little 
phase  change  to  drastically  alter  the  pattern.  The  effect  of  the 
out  of  nhase  element  is  to  focus  more  of  the  energy  in  the  lit  region 
of  the  cylinder.  By  changing  elements  2 or  3 such  that  they  are  out 
of  phase  with  each  other  should  eliminate  the  null  at  e = 0 and  also 
allow  the  lobes  to  be  moved  or  combined  into  a narrower  beam  as  well 
as  move  the  nulls. 
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Pig  10.  Antenna  pattern  of  monopole  on  ma.ior  axis  of 
elliptical  cylinder  (calculated  and 
experimental ) 


Fig.  11.  Roll  plane  pattern  (E.f)  for  a 1/25  scale  model 
of  a KC— 135  with  a a/4  monopole  on  the  fuselage 
forward  of  the  wings,  (from  Ref  14:64) 
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7.  Conclusions 


The  geometrical  theory  of  diffraction  has  been  extended  to 
treat  mono poles  on  perfectly  conducting,  cylinders  of  variable 
curvature.  Following  Pathak's  approach  (Ref  23:19-20)  of  defining 
a launch  coefficient  for  a slot  in  a cylinder,  a launch  coefficient 
for  a monopole  was  obtained.  For  the  monopole  case,  the  launch 
coefficients  are  the  same  as  for  an  axially  mounted  and  a length 
wise  mounted  slot  superimposed  on  each  other. 

By  carefully  keeping  track  of  the  distance  traveled  by  the 
diffracted  rays  an  interference  pattern  for  two  or  more  rays  or 
two  or  more  sources  was  obtained.  If  the  field  pattern  is  analyzed 
in  such  a way  as  to  lot  the  observation  distance  approach  infinity, 
then  the  phase  could  be  approximated  as  the  phase  at  the  diffraction 
point  for  the  ray  leaving  the  cylinder  ( ) • In  the  illuminated 
region  the  phase  would  be  assumed  to  be  that  of  the  source.  Although 
this  approximation  is  valid  it  was  assumed  for  this  problem  that 
since  this  antenna  is  to  serve  as  an  air  to  air  communications  link 
that  the  range  was  not  infinite.  A nominal  distance  between  com- 
munication points  of  300  meters  was  assumed  such  as  to  account  for 
aircraft  in  close  proximity  to  each  other.  The  field  further  out 
will  ju3t  have  a smaller  ripple  on  the  lobe  boundary  due  to  the  phases 
of  the  wavefronts  being  nearer  the  same  as  a function  of  their  dis- 
tance allowing  the  previous,  approximation  to  be  valid. 

The  effect  of  the  cylinder  is  to  broaden  the  beamwidth  in  the 
principle  plane.  By  mounting  the  ar.tenna  perrendi cular  to  the  sur- 
face, the  condition  of  the  antennas  no?  being  parallel  in  the  x,  y 
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plane  resulted.  Linear  array  theory  was  therefore  not  valid  (i.e., 
an  array  factor  is  not  present).  If  the  antenna  were  mounted  par- 
allel to  each  other,  the  effect  would  be  similar  (in  the  lit  region) 
to  tilting  a monopole  on  a flat  ground  plane  so  that  one  side  of  the 
beam  is  directed  into  instead  of  along  the  ground  plane.  Again, 
linear  antenna  theory  would  not  apply.  If  the  field  point  is  in  the 
deep  shadow  region  of  all  three  elements  of  the  array  then  the 
result  is  an  equivalent  linear  array  summation  at  the  diffraction 
point  Qg  but  it  is  only  valid  for  that  field  point  producing  a 
tangent  at  Qg.  It  can  also  be  assumed  that  the  angle  with  which  the 
antenna  is  attached  to  the  surface  and  not  its  angle  with  respect 
to  the  center  of  the  cylinder  controls  the  location  of  the  nulls. 

The  array  analyzed  provides  a broad  beam  for  communication 
problems  but  with  high  mobility  of  nulls  for  a null  steering  use. 

Due  to  the  creeping  wave  effect,  it  appears  it  would  be  very  hard 
to  create  a single  beam  in  one  general  direction.  The  symmetrical 
cylinder  would  tend  to  cause  multiple  lobes  as  in  Pig.  16.  An  end 
fire  or  broadside  pattern  in  the  x,  y plane  is  not  readily  ob- 
tainable because  the  elements  are  not  parallel.  However  when 
analyzing  the  fields  in  the  x,  z plane  (rt  = 0,  180)  the  antenna 
elements  appear  to  be  a. linear  array  with  antenna  2 and  3 reduced 

in  height.  Planes  of  view  in  between  (0  < « < 90)  would  be  a 

s 

combination  of  the  two  approachs  although  the  GTD  can  still  be 
applied  directly  and  obtain  valid  results. 
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Appendix  A 


Derivation  of  D and  «< 
m m 


The  first  step  is  to  formulate  and  solve  the  boundary  value 
problem  of  diffraction  of  a cylindrical  wave  by  a circular  cylinder 
of  radius  a.  Then  the  solution  will  be  expanded  asymptotically  for 
large  values  of  ka. 

In  order  to  simplify  computations  and  avoid  repetition,  a 
scalar  field  is  assumed  with  conditions  stated  to  allow  for  hard 
and  soft  boundary  conditions.  This  scalar  field  u(r,©)  must  satisfy 
the  equation 


(V2  + k2)u«$(r-/9  )i(«) 


(54) 


Two  different  boundary  conditions  at  r = a will  be  considered.  They 
correspond  to  different  physical  problems. 

By  expanding  the  diffracted  field  in  cylindrical  wave  functions 
and  imposing  the  boundary  conditions  ( Ref  24:Ch.  5)  the  field  for 
the  ray  OQ^C^P  as  shown  in  Fig.  17  becomes 


OJ  (ka) 

m -Hi2)(kr<) 


GH^^ka) 


f/2)(kr>) 
r > a 


(55) 


Here  r>  and 
ouantities 


r^  are, 

r and 


respectively,  the  larger  and  smaller  of  the 
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Pig.  17.  Geometry  for  Diffraction 

The  operator  G,  which  is  determined  by  the  boundary  conditions 
as  given  in  Table  II. 

Table  II 

The  Operator  £ 


G Boundary  condition  Field  Equivalent 


Case  I:  1 

u = 0 

s 

u = E ; perfectly 

conducting 

cylinder 

Case  II: 

in 

IT  = 0 
an 

h 

u = H ; perfectly 
z 

conducting 

cylinder 

The  representation  of  the  field  as  given  by  Eq.  (55)  is  valid 
only  for  small  values  of  ka.  A numerically  useful  result  for  large 
values  of  ka  is  obtained  by  emnloying  the  well-known  Watson  trans- 
formation which  leads  to  an  asymptotic  expansion  of  the  field. 
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Accessible  descriptions  of  the  method  are  given  in  References  25  and 

26. 


i. 


9 

Letting  C be  a contour  encircling  the  real  axis  as  shown  in 
Fig.  18,  the  sum  in  Eq.  (55)  becomes 


u = _ l(  exn(jff9  - jr)2  *4  kr>  ) 

Vc  •W”)  (ka) 

(56) 

• (GH^ka)^^)  - GJ^(ka)H^2)(kr<))djr 


When  yis  replaced  by  -Von  that  part  of  the  contour  for  which 
^ 0,  Eq.  (56)  becomes 

u.-i^ 


4 | D ain{VTl) 


GH^ka) 


• (Gf4r2)(ka)Jv(kr<)  - GJK(ka)H^2)(kr<))d»^  (57) 


V PLANE 

* 

f 

t 

% 

\ 

\ 

V 

« 

1 

* 

?itr.  18.  The  noth  of  inf erration  for  the 
Watson  type  integral  re  "aoentation  of  the 
field  diffracted  V cvlir.der 
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If  the  point  (r,©)  lies  in  the  geometrical  shadow,  the  contour 
D can  be  closed  in  the  lower  half-plane.  The  integral  can  now  be 
evaluated  by  computing  the  residues  at  the  zeroes  y of  GhL2  ka  )=0. 
This  yields 


■ . -f-: 


— (w.c  - "» Gii)(ka> 


m 

sin(VTT) 

ID 


7^H^2)(ka)  -m 

• m 


fl^.2)(kr)fi^2)(k/0;  (58) 


The  values  for  are  given  in  Table  IV,  Appendix  D.  Asymptotically 
expanding  the  Hankel  function  appearing  in  Eq.  (58)  (Ref  18:365-866): 


b' 

'-7:: 


U 


.1  «p(— .ikUr2  - a2)1/2  + (f>2  - a2)l/2)) 

2k(r^  - 77] 


cos  (^(e  - rr) 


’ exp(  cos-1  (-|)  + cos~x(|) ) ) — m,  n ' 

9 f-GH(/j(ka) 


/ 3-  > 


GH^^ka) 


*v  K 


m 


(59) 


Using  Euler's  equation  to  replace  cos(^(^  ~ ) ) yields 


s7n7 v~eY 
m 


. .1  e*p(-.‘k((r2  - a2)4/2  + (?2  - ag)l/2)) 

2k(/-a2)i',Y-a2)1''4 

. ^exp(-jtfri)  + exp(-,j]/my2) 
mTl  1 - exp(-j2n^) 

GH^Va) 

• exn(  jiT(cos-1(-p)  + cos_1(-|)))r  ^ 


W2W) 


(60) 
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where  Vj  3 6 and  ^2  = -£  Eq.  (60)  contains  the  leading  terms 

in  the  asymptotic  expansion  of  the  field  in  the  shadow  region. 

To  calculate  the  GTD  field,  Eq.  (12)  in  the  text  of  the  thesis 

* 

is  used  as  a starting  point.  For  a torsionless  surface  path  on  a 
cylinder,  no  E field  exists  in  the  £ direction  so  Eq.  (10)  is  equal 
to  zero  and  Eq.(l2)  may  be  written  as 


./dw  P 

u - * s> 


w ’ /®2  . 
m-1  /Qj  J 


\ - ■^=^exp(-j(kR  + «p) 
1 l/8kfi 


Due  to  the  symmetry  of  the  problem  it  is  clear  that  the  geo- 
desics are  arcs  of  the  generating  circle  of  the  cylinder.  Further- 
more, the  quantity  d^/di^  has  the  value  unity  since  a band  of  dif- 
fracted rays  does  not  spread  on  the  cylinder.  Also,  the  radius  of 

l/P 

curvature  p in  Eq.  (6l)  is  infinite  so  that  (yO/(y©  + s)s)  '*  assumes 

-l/2 

its  limiting  value  s ' . Moreover,  because  the  geodesic  curvature 

of  a ray  is  constant,  the  quantities  and  D^(  ) 3-i*e  constants 

with  Dm(Q^)  = ^(^2)*  introducing  these  simplifications  into  Eq.  (6l) 

2 2l/2 

and  using  Eq.  (62),  noting  that  3 at  is  (p  - a ) ' . Then 
Eq.  (6l)  becomes 


* A 

A torsionless  surface  path  is  one  where  the  direction  of  b 
does  not  change  along  a geodesic  rath  on  the  cylinder. 
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~^I^~ — exp(-j(k((p2  - a2)1/2  + s)  + <f)) 

- a2)1/2 

• ^D2(Q)exp(-jk  - (63) 

m=l 


Eq.  (63)  gives  the  field  associated  with  any  ray  from  0 to  P 
having  an  arc  length  X on  the  cylinder.  For  the  ray  OQ^QgF,  the 
length  X for  a circle  has  the  value  XQ  = a&  - a cos  ^(a/p)  - 


a cos~^(a/r),  while  s = (r^  - a^)L/c.  In  addition  to  this  ray, 

all  of  the  rays  which  are  tangent  to  the  cylinder  at  and  encircle 

the  cylinder  g times  before  emerging  at  Q2  are  also  diffracted  rays 

through  P.  For  these  rays  j has  the  value  X = / + 2gUa.  Inserting 

g o 

this  value  of  X into  Eq.  (63)  and  summing  over  g,  the  diffracted 
field  at  P due  to  those  rays  which  are  tangent  at  is  given  by 


2,1/2 


u = 


inf  - tPwtttk  - >2)1/2  * <r2  - ^4 

«o 

• ^D2exp((-jk  - <*  )/)(l  - exp(2ir(-jka  - a<X  J))-1  (64) 

. hi  m o ui 


Along  with  the  family  of  rays  described  above,  there  exists 
. other  rays  which  are  tangent  to  the  cylinder  at  Q^,  encircle  the 

cylinder  g times  and  then  emerge  tangentially  at  and  pass  through 

P.  Computing  the  field  along  these  rays  yields  a result  of  the  form 
of  Sq.  (64)  with  ©replaced  by  2TT-  ©in  Adding  this  result  to 

Eq.  (65),  we  obtain  for  the  total  diffracted  field  at  P 
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u 


-.1I^exp(-.jk UP2  - a2)1/2  + (r2  - a2)l/2)  + <f) 

ikif2  - a2)f/4(r2  - a2)1/4 


D^eip((-jka  - a*^)  + exp((-jka  - z\)Y2) 
1-  exp(2n(-jka  - a°<  ) ) 


m=l 


m 


17 


exp((jka  + a«m)(cos“1(-|)  + cos-1(jf))) 


(65) 


where  *i-*  and  ^2  s 2TT  - ©#  Comparing  Eq.  (65)  to  Eq.  (60)  and 
ignoring  the  amplitude  and  the  phase  of  the  source  gives 


-jk  + j&ra 


-1 


D2  = 
m 


-ft) 


-1/2  GH^ka) 

1 ' in 


(66) 


Expansion  of  GH^  ^(ka^/f^/JjfjGH^ 2 ^(ka  jj  is  done  in  Appendix  D. 


m "m 

Voltmer's  generalizations  for  surfaces  of  variable  curvature 


are  contained  in  Table  I in  the  text,  and  are  seen  to  be  dependent 
on  not  only  a . the  radius  of  curvature  of  the  surface  with  respect 

o 

to  arc  length  along  the  ray  trajectory,  hut  also  and  /^n» 

where  the  dot  denotes  a derivative  with  respect  to  arc  length  and 

is  the  radius  of  curvature  of  the  surface  in  the  direction  of  the 
binormal  to  the  ray. 


■ 


r 


- 1 >» 


Appendix  B 

The  Launching  Coefficient 

The  launching  coefficinet  is  derived  by  asymptotically  expan- 
ding a canonical  problem  for  ka  large  and  comparing  this  to  the  GTD 
solution  as  was  done  for  the  diffraction  coefficients  in  Appendix  A. 
Start  with  Eq.  (58)  in  Appendix  A. 


co 


,(l)f 


Cos(v  (®  “ ir))  GV  (o)  (o) 

■ - fS  A,, 


(58) 


iv  v , 


m 


cos(ym(©  - n)  exp(-jj/m'f/1)  + exp(-jv^^) 

3eplacing  ' s~nT7  rf) = 1 exp(-7*V ’ and  letting 


(l  - exp(-j27T^r ) )“^  a 1 for  ka  large  gives 

,(1) 


..-■yZ  -/vf^)-i2)(^)i2) 

4 mTl  GH  3 (ka)  ^m  *m 


expi-j^)  + exp(-jymy2) 


(y) 


(67) 


where  ^ = © and  V^>  = 2 TT  — d.  Using  the  definition  of  the  diffrac- 
tion coefficient  from  Arnendix  A 


(68) 


*42)(k/»)  I^xp(v  ■ (2/+ 1)5) 


and  defining  the  boundary  conditions  as 


Y>  = 0 for  p = a (hard  boundary)  (69) 


and 


(soft  boundary)  (70) 


u(a^S^)  = 0 

then  for  the  hard  boundary  condition,  noting  the  results  for  a point 
source  in  Eq.  (23) 

a 


- Ix‘fu(?»xa) 


Then  Eq.  (67)  becomes 


jkl 


idn k j m=l 


V 

m 


(71) 


[exp(-j^)  + e:,p(-jyj2)] 
~h  expf-.jk.^ 

* Tp' 


(72) 


Heferring  to  Fig.  17*  the  GTD  solution  for  E with  a launching 
coefficient  replacing  the  field  incident  at  is 


ra=l 


»xp((-:ka  - Yp'  * 


(73) 


i 
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Comparing  Eqs.  (72)  and  (73)  gives 


°C  = ka  " (?5) 


Replacing  the  hard  coefficients  with  the  soft  coefficients  and 

(2) 

differentiating  '(ka)  with  respect  to  ka  will  solve  the  problem 

m (2)  (2) 

for  the  soft  boundary  condition.  Hv  '(ka)  and  Hv  ''(ka)  are  ex- 

Ki  ^m 

panded  in  Appendix  D. 

Comparing  this  to  Pathak's  work  on  slots  in  cylinders  (Ref  23: 
46-62)  shows  that  the  same  launch  coefficient  results  for  a magnetic 
current  source  in  the  z direction.  The  soft  boundary  condition  for 
a monopole  is  equivalent  to  a launch  coefficient  of  a magnetic  line 

Jk  V 

source  oriented  in  the  ©direction  (circumferential  slot).  The 
monopole  can  then  be  modeled  as  2 magnetic  line  sources  perpendicular 
to  the  direction  of  the  monopole  electric  source  for  fields  in  the 
deep  shadow  region. 


Appendix  C 


Airjr  Functions,  Derivatives  and  Zeros' 


The  Airy  differential  equation, 


d w(z)  / \ 

t - zw(  z ) = 0 


(76) 


dz 


yields  two  independent  solutions  which  assumes  several  forms. 
The  form  used  in  tbi3  thesis  is  defined  by  Miller  as 


Ai(z)  = j exp(zt  - t3/3)dt 


(77) 


s 


which  for  z =■  x (x  real)  may  be  written  as 


Ai(x)  = = ( cos(xt  + t^/3)dt 

jo 


(78) 


is  shown  in  Figure  19.  A second  independent  solution  is  given  by 


-j2H/3 


Ai(ze  J 


if  . 


2WjJre::p(ze 


t - tV3)dt 


,.1211/3  / 3 

-2Tfj — I exp(  zt  - tJ/3  )dt 

r n 


(79) 


while  the  l-.ronskian  relation  satisfied  by  these  solutions  is 
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J 


(80) 


2 n 


Ai(z)Ai  '(ze~^21T/^)e-^2lT^  - Ai( ze“^2TT/^)Ai  ' ( z) 


The  derivatives  are  obtained  by  differentiating  under  the  integrals 
of  Eqs.  (77)  and  (79). 

The  roots  of  the  Airy  function  and  its  derivative  which 
satisfy 

Ai(z)  - 0 and  Ai'(z)  * 0 

are  real  and  negative  and  are  denoted  as  -q^  and  -q^,  respectively, 
where  m denotes  the  number  of  the  root.  The  first  two  roots  are 
presented  in  Table  III  as  are  the  values  of  Ai(-q"m)  and  Ai'(-qm). 

Table.  III.. 

Airy  Functions,  Derivatives,  and  Zeros  • 


m 


m 


1 +2.33311 

2 +4.08795 


m 


#i’ 


Ai<-V 


+1.01879 

+3.24820 


+0.70121  +0.53566 

-0.80311  -0.41902 

(from  Ref  20:149) 
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Appendix  D 


Hankel  Functions 


An  evaluation  of  the  integral  resulting  from  the  Watson  trans- 
formation requires  the  calculation  of  the  roots  of  the  equation 


r<2) 


Gil;,  '(ka)  = 0 


(81) 


where 


Soft  EM  cylinder 


Hard  EM  cylinder 


JX  IX  = kP 

A closed-form  evaluation  of  these  zeros  is  not  possible.  However, 
an  asymptotic  representation  of  the  roots  may  be  obtained  for  large 
values  of  ka.  The  procedure  is  outlined  in  Reference  20,  pages 
127  to  133.  The  results  are  listed  in  Tables  IV  and  V on  the  next 
page.  The  asymntotic  approximations  for  the  Hankel  functions  of  the 
second  order  are 


1/3 


1 uc-s) 


i . 


• 2T7 

rT 


,ka/ 


i - (ff1 


,2V  \ 

* 


M/  15 


(82) 

(83) 


an  alternative  form  is 


r'W  -s(kP-  (2K*  lj4! 


(65) 
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Table  IV 


Hankel  Function  Zeros, 


GH^2)(ka) 


0 


Surface 


Soft 


Hard 

Cylinder 


Zeros- V 


(from  Ref  20;145) 


Table  V 


GH^(ka) 

Hankel  Function  Quotient,  r tt?\ 

bi(ka) 


Surface 


G-Crerator 


Quotient 


Soft 


Hard 


l 

d(*0) 


e~j5iy6  2&  (1  (lX/l  S e-d*/3 

2*T( Ai  • (-q  ))2  U W 30 


m 


.-^/6(f)Vi  .2/3  £ /3 

r n2  (ka/  30  e 


(from  Ref  20:136) 
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Appendix  E 


The  Field  in  the  Transition  Region 


The  material  for  this  Appendix  is  taken  in  part  from  Reference 
23,  pages  63  to  70.  In  the  transition  region  adjacent  to  the 
shadow  boundary,  the  representation  of  the  field  radiated  by  the 
monopole  described  by  the  launching  coefficient  in  Appendix  B is 
poorly  convergent  so  a different  form  is  needed.  It  will  be  seen 
that  a suitable  asymptotic  solution  can  be  obtained  in  terms  of 
integrals.  The  first  order  approximation  is  adequate  for  ka  > 3, 
and  since  the  added  complexity  is  hard  to  justify,  only  the  first 
order  will  approximation  will  be  derived. 

For  the  principle  plane  analysis  of  the  cylinder  (looking 
length-wise  down  the  cylinder),  only  the  hard  boundary  condition 
exists,  therefore  only  it  will  be  treated.  The  soft  condition 
is  solved  identically  as  below  except  without  taking  the  iu/j^. 

The  Watson  transform  is  a~ain  arolied  to  Eq.  (55),  but  this 
time  the  resulting  expression  is  left  in  the  form  of  an  integral. 

In  the  far  zone,  H^,  '(kp)  is  approximated  asymptotically  as  in  Eq. 
(68)  so  that  (see  Ref  18:865-366) 


^ u -k  exp(-jk^)^  [ djA 

V 2f8nk  j*  TiT  )_  oo 

^'(kajH^ka)  - ^^(kajH^ka) 
H(^(ka)  ” 

. exp(-jT/(V|  + 2W|))  + exp(-jV(V.>  + 2TT-f) 


?rom  the  Wronskian  relationship  for  the  Hankel  function. 


Hi1,,<ka»Ii2,<ka>  - Hi2)’(ka)H!1,(ka)  -in k 


exp(-jk/ 


*o  />' 

K 


bp  kaV8nkj’  Vf~'  ^ )-  ao 

(exp(-jm(V>1  + 2T ti))  + exp(-jm(^  + 2ttf)) 

42W) 

Replacing  K^(ka)  by  its  Watson  approximations  (for  m*ka»l) 

m 

gives  rise  to  the  Fock  function  representations  for  the^u/^?: 

00  r 2/3 

T(  ^0- » (i? > ’ 

f exp(-j£jT) 

(-3 1 In)  yp  — 'w^F7 dr  + exp(-3ka(f2  + 2JTf)) 

■ kaY3,  ,._Y  exp(-J<2r)  I.*,,.*,, 

<T>  S^r! — 4TJ  ^ (l 

.here  H<2)(ka)  ~^§j)  »2W*  and  ^ ' (¥> 

Here  ’.^(T)  is  the  Fock-type  Airy  function  related  to  the  Killer- 
type  airy  function  by 

w2(T)  = 2Vrexp(-jTT/6)Ai(-Texp(  jrf/3))  (, 

or  = ||,  |»exp(rz  - z3^3)dz  ^ 

with  the  contour  of  integration  shown  in  Appendix  C . 

Introducing  the  Fock  functions 
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then  Eq.  (87)  becomes 


>u 


g^-^expt-jka^)  + g($2)  . 

• expf-okafc) 

If 


(91) 


Where  only  the  f=Q  term  is  retained  for  large  ka.  Substituting 
this  into  Eq.  (71 )»  the  far-zone  electric  fields  for  the  line  source 


^ T • 
x I is 

X 


E 'v-  C I 

© x 


K>  /• 

. (gt^1)exp(-jkaY/1)  + g(^2)exp(-3ka'/'2) ) . 


exp  (-JM) 

' *7 

Comparing  with  Eq.  (73  ) yields 


(92) 


▼(Q1*Q2)  - g(^1)exp(-jka'^)exp(-jk^) 

Tabulation  of  g(^)  is  included  in  Appendix  G. 


(93) 


distance  of  300  meters.  The  angle  v is  the  angle  the  field  radius 
vector  makes  with  the  verticle — the  same  as  the  parametric  v used 
in  the  mathematical  model  of  the  cylinder.  The  field  is  calculated 
for v between  0 and  l80°  in  steps  of  5°«  Due  to  symmetry  the  field 
due  to  antenna  2 in  the  180°  to  360°  range  is  the  same  as  for 
antenna  3 in  the  0 to  180°  range  with  the  direction  of  the  field 
vector  turned  around  (shifted  180°  spatially). 

To  calculate  the  total  fields  using  these  tables,  add  the 
&■-  source  phase  of  the  antenna  in  question  to  that  in  the  chart  for 

that  antenna,  and  multiply  the  x and  y quantities  by  the  source 
amplitude.  Sow  for  each  field  point  determined  by  v,  an  x,y, phase 
data  point  for  each  antenna  is  taken  from  the  table.  To  obtain  the 
total  field  these  three  vectors  are  added  vectorially. 

Also  noted  are  the  equations  which  were  used  to  calculate  the 
values  in  the  chart.  The  equations  are 

Eq.  # Equations  used 

1 Lit  region  Eq.  (40) 

2 Lit  side  of  transition  region 

' Eq.  (41) 

3 Shadow  side  of  transition  region 
Sq. (42) 

4 Deep  shadow  region  Eq.  (43) 
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Table  VI 


Field  Data 


Antenna  1 


V 

Eq.  # 

Magnitude 

X 

y 

Phase J-j(  • 

0 

1 

0.000 

0.000 

0.000 

0.00 

5 

1 

2.944 

-0.258 

2.933 

2.76 

10 

1 

5.865 

-1.024 

5.766 

11.01 

15 

1 

8.740 

-2.276 

8.438 

24.70 

20 

1 

11.545 

-3.957 

10.798 

43.70 

25 

1 

14.259 

—6 . 06  2 

12.906 

67.89 

30 

1 

16.861 

-8.479 

14.574 

97.07 

35 

1 

19.330 

-11.147 

15.791 

131.00 

40 

1 

21.646 

-13.986 

16.521 

169.44 

45 

1 

23.793 

-16.903 

16.744 

212.09 

50 

2 

24.000 

-18.463 

15.337 

258.61 

55 

2 

24.510 

-22.623 

15.659 

306.28 

60 

2 

25.068 

-21.781 

12.409 

3^4.02 

65 

2 

25.479 

-23.166 

10.627 

28.17 

70 

2 

25.555 

-23.907 

8.589 

76.30 

75 

2 

25.480 

-24.654 

6.436 

127.84 

80 

2 

25.026 

-24.674 

4.182 

182.32 

85 

2 

24.304 

-22.965 

1.958 

239.15 

90 

23.316 

-23.316 

0.000 

3.06 

95 

3 

22.223 

-22.138 

-1.937 

79.95 

100 

3 

21.082 

-20.762 

-3.661 

156.54 

105 

3 

19.826 

-19.150 

-5.131 

231.73 

110 

3 

18.771 

-17.639 

-6.420 

303.32 

115 

3 

17.618 

-15.967 

-7.^46 

3.05 

120 

3 

16.366 

-14.173 

-3.163 

73.22 

125  ’ 

3 

15.310 

-12.541 

-3.781 

137.67 

130 

4 

14.466 

-11.082 

-9.299 

195.93 

135 

4 

11.399 

-8.060 

-8.060 

208.36 

140 

4 

11.468 

-7.372 

-8.785 

284.12 

145 

4 

10.830 

-6.212 

-8.871 

344.67 

150 

4 

7.950 

-3.975 

-6.885 

50.41 

155 

4 

9.520 

-4.023 

-8.628 

107.02 

160 

4 

6.316 

-2.160 

-5.935 

128.98 

165 

4 

8.530 

-2.221 

-8.238 

164.79 

170 

4 

3.831 

-O.665 

-3.773 

145.14 

175 

4 

7.329 

-0.632 

-7.800 

164.79 

180 

4 

0.000 

0.000 

0.000 

0.00 

... 1 


Table  VII 


Field  Data 
Antenna  2 


V 

Eq.  # 

Magnitude 

X 

7 

Phase[-j(«  jj 

0 

1 

15.418 

0.000 

15.418 

22.90 

5 

1 

17.970 

-1.591 

17.899 

33.46 

10 

1 

20.382 

-3.578 

20.066 

49.28 

15 

1 

22.635 

-5.912 

21.849 

70.24 

20 

1 

24.712 

-8.521 

23.196 

96.17 

25 

1 

25.595 

-10.897 

23.159 

126.87 

30 

2 

25.222 

-12.699 

21.792 

156.53 

35 

2 

25.220 

-14.555 

20.596 

186.01 

40 

2 

25.546 

-16. 5 16 

19.48 9 

220.88 

45 

2 

25.606 

-18.201 

18.012 

260.81 

50 

2 

30.279 

-19.442 

16.135 

305.37 

55 

2 

24.678 

-20.296 

14.038 

344.71 

60 

2 

23.830 

-20.707 

11.792 

36.77 

65 

3 

22.866 

-20.724 

9.664 

86.77 

70 

3 

21.723 

-20.413 

7.430 

132.00 

75 

3 

20.575 

-19.374 

5.325 

234.78 

80 

3 

19.410 

-19.115 

3.371 

7.94 

35 

3 

18.229 

-18.160 

1.589 

83.45 

90 

3 

17.083 

-17.083 

0.000 

196.85 

95 

3 

15.960 

-15.899 

-1.391 

274.39 

100 

4 

14.854 

-14.628 

-2.579 

352.03 

105 

4 

10.921 

-10.590 

-2.827 

296.08 

110 

4 

10.639 

-9.997 

-3.639 

18.86 

115  ' 

4 

10.469 

-9.483 

-4.424 

110.64 

120 

4 

9.359 

-8.538 

-4.930 

190.99 

125 

4 

8.484 

-6.950 

-4.866 

272.01 

130 

4 

8.274 

-6.338 

-5.318 

347.50 

135 

4 

6.569 

-4.645 

-4.645 

79.73 

140 

4 

7.771 

-4.995 

-5.953 

144.93 

145 

4 

6.136 

-3.519 

-5.026 

189.05 

150 

4 

3.825 

-1.912 

-3.313 

286.93 

155 

4 

7.305 

-3.037 

-6.621 

315.70 

160 

4 

1.370 

-0.469 

-1.237 

325.07 

165 

4 

7.241 

-1.374 

-6.994 

24.83 

170 

4 

6.601 

-1.146 

-6.501 

59.56 

175 

4 

8.130 

-0.709 

-8.099 

9.34 

180 

4 

5.024 

0.000 

-5.025 

117.38 

j 
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Table  VIII 


Field  Data 
Antenna  3 


V 

Eq.  # 

Magnitude 

z 

y 

Phase  [- j (•  )J 

0 

1 

15.418 

0.000 

-15.418 

22.90 

5 

1 

12.745 

1.110 

-12.697 

17.67 

10 

1 

9.973 

1.732 

-9.821 

17.83 

15 

1 

7.124 

1.844 

-6  • 88l 

23.36 

20 

1 

4.440 

1.443 

-3.966 

34.23 

25 

1 

1.285 

0.543 

-1.165 

50.35 

30 

1 

1.658 

-0.829 

1.436 

71.60 

35 

1 

4.587 

-2.631 

3.758 

97.81 

40 

1 

7.478 

-4.827 

5.712 

128.78 

45 

1 

10.308 

-7.319 

7.260 

164.27 

50 

1 

13.056 

-10.039 

8.349 

204.00 

55 

1 

15.700 

-12.905 

8.943 

247.67 

60 

1 

18.220 

-15.827 

9.027 

294.94 

65 

1 

20.596 

-18.715 

8.600 

345.45 

70 

1 

20.810 

-21.480 

7.678 

38.80 

75 

1 

24.846 

-24.039 

6.283 

94.60 

80 

1 

25.670 

-25.313 

4.300 

152.40 

85 

2 

25.797 

-25.713 

2.082 

210.15 

90 

2 

25.842 

-25.841 

0.000 

268.30 

95 

2 

25.972 

-25.858 

-2.431 

327.97 

100 

2 

26.147 

-25.719 

-4.707 

11.99 

105 

2 

25.405 

-24.496 

-6.734 

39.95 

110 

2 

24.700 

-23.157 

-8.594 

150.92 

115 

2 

22.770 

-20.577 

-9.749 

212.65 

120 

3 

21.622 

-18.725 

-10.811 

289.60 

125 

3 

21.618 

-17.708 

-12.400 

302.50 

130 

3 

20.368 

-15.603 

-13.092 

0.72 

135 

3 

19.219 

-13.590 

-13.590 

58.52 

140 

3 

18.288 

-11.755 

-14.009 

103.95 

145 

3 

16.926 

-9.708 

-13.865 

161.14 

150 

3 

15.697 

-7.848 

-13.594 

206.92 

155 

4 

14.397 

-6.084 

-13.048 

273.65 

160 

4 

11.630 

-3.978 

-IO.929 

304.16 

165 

4 

11.513 

-2.981 

-11.126 

347.31 

170 

4 

3.274 

-1.437 

-8.148 

329.34 

175 

4 

10.593 

-0.923 

-10.553 

345-40 

180 

4 

3.329 

0.000 

-3.291 

290.14 
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Appendix  G 


1 


Table  of  the  Hard  Pock  Function 


The  function  is  defined  as 


g(C) 


exp(-j(t> 


w{(t) 


dt 


Table  IX 

The  Hard  Fock  Function 


6 

Modulus 

Argument 

-1.00 

1.861 

15-460 

-0.95 

1.847 

12.626 

-0.90 

1.833 

10.077 

-0.85 

1.818 

7.801 

-0.80 

1.802 

5.785 

-0.75 

1.785 

4.016 

-0.70 

1.766 

2.482 

-0.65 

1.747 

1.174 

-0.60 

1.726 

0.079 

-0.55 

1.705 

359.187 

-0.50 

1.682 

358.487 

-0.45 

1.658 

357.971 

-0.40 

1.633 

357.627 

-0.35 

1.607 

357.446 

-0.30 

1.580 

357.419 

-0.25 

1.552 

357.536 

-0.20 

1.523 

357.790 

-0.15 

1.493 

358.172 

-0.10 

1.463 

358.672 

-0.05 

1.431 

359.284 

0.00 

1.399 

0.000 

0.05 

1.367 

0.812 

0.10 

1.334 

1.713 

0.15 

1.300 

2.6  96 

0.20 

1.266 

3.755 

0.25 

1.232 

4.883 

0.30 

1.197 

6.074 

0.35 

1.163 

7.323 

0.40 

1.128 

3.624 

0.45 

1.093 

9.071 

0.50 

1.059 

11.360 

0.55 

1.025 

12.786 

65 

I 


0.60 0.991 14.244 
0.65  0.958  15.731 
0.70  0.924  17.243 
0.75  0.892  18.776 
0.80  0.860  20.326 
0.85  0.828  21.893 
0.90  0.798  23.474 
0.95  0.768  25.070 
1.00  0.738  26.683 


(from  Ref  21sch.  4) 
For  more  accurate  values  see  Ref  21. 
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